Cosmology with interaction in the dark sector 
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Unless some unknown symmetry in Nature prevents or suppresses a non-minimal coupling in the 
dark sector, the dark energy field may interact with the pressureless component of dark matter. In 
this paper, we investigate some cosmological consequences of a general model of interacting dark 
matter-dark energy characterized by a dimensionless parameter e. We derive a coupled scalar field 
version for this general class of scenarios and carry out a joint statistical analysis involving SNe la 
data (Legacy and Constitution sets), measurements of baryon acoustic oscillation peak at z = 0.20 
(2dFGRS) and z = 0.35 (SDSS), and measurements of the Hubble evolution H{z). For the specific 
case of vacuum decay {w — —1), we find that, although physically forbidden, a transfer of energy 
from dark matter to dark energy is favored by the data. 

PACS numbers: 98.80.Cq; 95.36.-l-x 
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I. INTRODUCTION 

As is well known, there is mounting observational ev- 
idence that our Universe is presently dominated by two 
exotic forms of matter or energy. Cold, nonbaryonic 
dark matter (CDM), which accounts for ~ 30% of the 
critical mass density and whose leading particle candi- 
dates are the axions and the supersymmetric charginos 
and neutralinos, was originally proposed to explain the 
general behavior of galactic rotation curves that differ 
significantly from the one predicted by Newtonian me- 
chanics. Later on, it was also realized that the same 
concept is necessary for explaining the evolution of the 
observed structure in the Universe from density inhomo- 
geneities of the size detected by a number of Cosmic Mi- 
crowave Background (CMB) experiments. Dark energy 
or quintessence, which accounts for ~ 70% of the crit- 
ical mass density and whose leading candidates are a 
cosmological constant A and a relic scalar field $, has 
been inferred from a combination of astronomical obser- 
vations which includes distance measurements of type la 
supernovae (SNe la) indicating that the expansion of the 
Universe is speeding up, CMB anisotropy data suggest- 
ing — Ij,^'^^ clustering estimates providing Qm — 0.3 
(see, e.g., [l| for current reviews). 

Another interesting aspect related to these two dark 
components is that, unless some unknown symmetry in 
Nature prevents or suppresses a non-minimal coupling 
between them (see for a discussion), the dark matter 
and dark energy fields may interact between themselves, 
giving rise to the so-called models of coupled quintessence 
(see H, 0, H, [1] and Refs. therein). From the observa- 
tional side, no piece of evidence has been so far unam- 
biguously presented against such an interaction, so that a 
weak coupling still below detection cannot be completely 
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excluded. In the cosmological context, models of coupled 
quintessence are capable of explaining the current cosmic 
acceleration, as well as other recent observational results 
01 S @|- From the theoretical viewpoint, however, cri- 
tiques to these scenarios do exist and are mainly related 
to the fact that in order to establish a model and study 
its cosmological consequences, one needs first to specify 
a phenomenological coupling F between the cosmic com- 
ponents. 

In this paper, instead of adopting the traditional ap- 
proach, we follow the qualitative arguments used in 
Refs. 0, [^, and deduce the interacting law from a 
simple argument about the effect of the dark energy de- 
cay on the dark matter expansion rate. As a step towards 
fundamental physics motivation for this phenomenolog- 
ical dark matter/dark energy interaction, we study how 
this class of coupled cosmologies can be interpreted in 
terms of classical scalar field dynamics. We also discuss 
current observational constrains on the interacting term 
e. To this end we use two different samples of SNe la 
data, namely, the SNLS collaboration sample of 115 data 
points [ll| and the most recent SNe la compilation, the 
so-called Constitution set of 397 SNe la covering a red- 
shift range from z = 0.015 to z = 1.551 [l^. Given the 
complementarity of SNe la data with other cosmologi- 
cal observables, a joint analysis with measurements of 
baryon acoustic oscillation (BAG) peak at z = 0.20 [l3| 
and z = 0.35 11411 . and measurements of the Hubble evo- 
lution H{z) [l5| is also performed. 



II. INTERACTING DARK MATTER/DARK 
ENERGY 

If dark energy and dark matter interacts, the energy 
density of this latter component will dilute at a differ- 
ent rate compared to its standard evolution, pdm oc a~^, 
where a is the cosmological scale factor. Thus, the de- 
viation from the standard dilution may be characterized 
by a constant e, such that 
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FIG. 1: Left: The ratio pde/pdm as a function of the redshift parameter z [Eqs. ([1} and ([3])] for some selected values of e and 
w and r2tie,o/f2dm,o — 3. Note that, for large and positive values of e, the relative contribution of dark energy to dark matter 
may be significant at early times. Right: The potential U(^) [Eqs. (O-®] as a function of the field for some selected values 
of w and e = 0.1 and Q,dm = 0.23. In this Panel, "NB" means that the baryons participate of the interacting process whereas 
"WB" means that they were regarded as a separetely conserved component contributing with Q,b,o — 0.0416. 



where we have set the present-day value of the cosmo- 
logical scale factor ao = 1. By assuniing that the radi- 
ation and baryonic fluids are separately conserved, the 
energy conservation law for the two interacting compo- 
nents reads UaT^^ = 0, where T = -I- T^f or, 
equivalently (in our notation, pm = pb + Pdm), 
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By combining the above equations and considering that 
the dark energy component is described by an equation 
of state parameter w = pdel Pde, one can show that the 
energy density of the dark energy component evolves as 
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where the integration constant pde,o is the present-day 
fraction of the dark energy density. Clearly, in the ab- 
sence of a coupling with the CDM component, i.e., e = 0, 
the conventional non-interacting quintessence scenario is 
fully recovered. For w; = — 1 and e 7^ 0, we may identify 
Pdefi = pv,o (the current value of the vacuum contribu- 
tion), and the above expression reduces to the vacuum 
decaying scenario recently discussed in Refs. 0, H, @- 

Fig. (la) shows the ratio Pde/pdm as a function of the 
redshift parameter [Eqs. ((T]) and ^] for some selected 
values of e and w and ^de,o/^dm.o — 3. It is shown that 
this fraction is a very sensitive function of the interacting 
parameter e. For instance, at 2; = 5, the dark energy con- 
tribution relative to pdm for e — 0.15 is about 75% larger 
than the same contribution for e — 0.10. This amounts to 
saying that for relatively large and positive values of e the 
energy density of the pressureless component is not neg- 
ligible at high- 2;. For the sake of comparison, the ACDM 
case (e = and w = —1) is also shown (solid line). 



III. SCALAR FIELD DESCRIPTION 

In order to find the scalar field counterpart for the flat 
interacting dark matter/energy model considered here 
(that shares the same dynamics and temperature evo- 
lution law), we use the procedure originally proposed in 
Ref. [1^ (see also Q for a general description with arbi- 
trary curvature). 

First, let us define the parameter 7* = —2H/3H^ = 
(1 -j-Loflde), which is just another way of writing the field 
equations 



87rG(pf, + Pdm + Pde) = SH^ 



inGpde = -2H - 3H^ , 



(4) 
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where pb, Pdm and pde are the energy densities of 
baryons, cold dark matter and dark energy, respectively, 
whereas pde stands for the dark energy pressure. Thus, 
any cosmology dynamically equivalent to the coupled 
quintessence scenario considered here must have the same 
7*. Following standard lines, we replace the energy den- 
sity and pressure in Eqs. (|4]) and (O by the correspond- 
ing scalar field expressions, i.e., pde P<p and pde —>■ P4>, 
where p^ — 0^/2 + V^(0) and p^ — (j)'^/2 - F(0) are, 
respectively, the energy deirsity and pressure associated 
with the coupled scalar field whose potential is 1^(0). 

By defining a new parameter x = 0^/(0^ -I- pm) with 
< X < 1, we can mairipulate the above equations to 
obtain (see for more details) 
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FIG. 2: Contours of x in planes w — e (left), f2dm,o — e (middle) and Odm.o — e with w = — 1 (right). These contours are 
drawn for Ax^ = 2.30 and 6.17. In all Panels, dashed lines correspond to the joint analysis involving SNLS + BAO + H(2) 
measurements whereas the solid ones to CS + BAO + H(z) data. The shadowed area in the Panel at right stands for the 
thermodynamical constraint on e discussed in Ref. [1]. 



which link directly the field and its potential with the 
related quantities of the coupled quintesence case. From 
Eq. (|6a|) , one can show that in terms of the scale factor 
a the field is given by 



da 
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where 7* (a) can be obtained from Eqs. H]) and ([3]), i.e., 



nb.oa-'^ + (1 + w)nde,oa-^^^+-''^ + Andm,oc 
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with A - ^il^ and B - 



(8) 

„i I toi^^ ^ — — 1 — . From now on, we 
consider x to be a constant, which is equivalent to im- 
pose the condition that the scalar field version mimics 
exactly the interaction rate of its phenomenological cou- 
pled counterpart [31 . From Eqs. ([Bb]) and ([5]), we also 
obtain 



U{a) 



-A - wB 
2 



-{1 + w) — w 
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where U — 8nGV{a)/3HQ and flde,o = 1 ^ ^b.o — 
3|w|f^dm,o/(3|w| - e). 

Figure (lb) shows the potential U{(j>) obtained from 
a numerical combination of Eqs. (l7])-(l9]) by assuming 
e = 0.1 and ^dm.o — 0.23. Four different cases are 
shown, namely, quintessence (w = —0.8) and phantom 
(w = —1.2) EoS with and without a separately conserved 
baryonic component flb^o- As one may also check, the po- 
tential shown in Fig. (1) reduces to 



24, 24> 
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when w — —1, where Ci, C2 and C3 are constants IQJ. In 
general, we can see that the effect of ~ 4.16% of baryons 
on the potential Uicj)), as well as a changing in the EoS 



values, produces only a small shift relative to the case in 
which the baryonic content participates of the interact- 
ing process or to the case of vacuum decay, so that the 
general potential for these scenarios belongs to the same 
class of potentials given analytically by Eq. PU)) [T7| . 



IV. CURRENT OBSERVATIONAL 
CONSTRAINTS 

The description of the interacting scenario discussed in 
Sections II and III clearly shows that it comprises a multi- 
tude of cosmological solutions for different combinations 
of e, w and ^dm,o- In this Section we investigate observa- 
tional bounds on the parametric spaces e—w and e—fldm.o 
from statistical analyses involving three classes of cosmo- 
logical observations. In our analysis we fix fli, — 0.0416 
from WMAP results [l^ , a value in good agreement with 
the constraints derived from primordial nucleosynthesis 

We use current distance measurements to SNe la, 
namely, the SNLS collaboration sample of 115 points [ll| 
and the most recent SNe la compilation, the so-called 
Constitution set (CS) The SNLS sample includes 

71 high-z SNe la in the redshift range 0.2 < z < 1 and 
44 low-z events compiled from the literature but ana- 
lyzed from the same manner as the high-z sample. The 
CS of 397 SNe la covers a redshift range from z = 0.015 
to z = 1.551, including 139 SNe la at z < 0.08, and con- 
stitutes the largest SNe la luminosity distance sample 
currently available. 

The BAO distance ratio from 2dFGRS and SDSS mea- 
surements is given by (0.35)/i:)y(0.2) = 1.812 ± 
0.060 [3, where 
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and H{z) is given by Eqs. (H]), ([3]) and ([1]). In our anal- 
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ysis, we also use this ratio along with 9 determinations 
of the Hubble evolution H{z), as given in Ref. [IB] (for 
more details on statistical analysis involving these data 
sets we refer the reader to Ref. [285). 

The results of our statistical analyses are shown in 
Fig. (2). As usual, the total likelihood is written as 
C oc where the function takes into account 

the three data sets discussed above. By marginalizing C 
over Qdm.o i'-^): we can quantify how much the plane e — uj 
(e — can be constrained by the data. The contour 
levels for these analyses are displayed, respectively, in 
Panels (2a) and (2b). In both Panels, dashed lines corre- 
spond to la and 2a values of the joint analysis involving 
SNLS + BAO + H(2;) measurements whereas the solid 
ones to CS + BAO + tl{z) data. Although the bounds on 
e from the latter analysis (CS) are tighter than those from 
the former one (SNLS), the important point worth em- 
phasizing here is that both negative and positive values 
for the interacting parameter are allowed by these anal- 
yses. Physically, this amounts to saying that not only an 
energy flow from dark energy to dark matter (e > 0) is 
observationally allowed but also a flow from dark matter 
to dark energy (e < 0) [see Eq. At 95.4% (C.L.), 

we have found w = -0.92 ± 0.15 and e = -OM^o'^l and 
w = -0.78 ± 0.08 and e = -IMtoH, respectively, for 
the SNLS -I- BAO + H(z) and CS + BAO + H(z) joint 
analyses. 

In Fig. (2c) we show a similar analysis by fixing the 
dark energy EoS at w = —1, which is fully equivalent to 
the vacuum decay scenario discussed in Refs. 0,i. As 
shown in Ref. [1], for the vacuum decay process the in- 
teracting parameter e is restricted by thermodynamical 
arguments to be positive [25|| . By considering this phys- 
ical constraint on e we have found at 95.4% (C.L.) e — 
O.Otom and f2d„i,o = 0.21 ± 0.06 [SNLS + BAO + H(z)] 
and e = O.OtJJ;^^ and rtdm,o = 0.24 ± 0.04 [CS + BAO -I- 
H(2;)], which is equivalent to the standard ACDM case. If 
we drop the thermodynamical constraint (e > 0), we find, 
respectively, e = -0.46to;e8 and il.dm,o = 0.19 ±0.05 and 
e = -0.24j:^ j° and Qdrnfi = 0.23 ±0.04 at 95.4% (C.L.). 



V. FINAL REMARKS 

We have discussed some cosmological consequences of 
an alternative mechanism of cosmic acceleration based on 
a general class of coupled quintessence scenarios whose 
interaction term is deduced from the effect of the dark 
energy on the CDM expansion rate. The resulting expres- 
sions for the model are parameterized by a small param- 
eter e and have many of the previous phenomenological 
approaches as a particular case. 

Following the procedure proposed in Refs. j^, [lB|, we 
have also discussed a scalar field description for these cos- 
mologies and found that this class of models is identified 
with a coupled quintessence field with the potential given 
byEqs. dZ!)-® [see also Fig. (1)]. From the observational 
point of view, we have investigated the current bounds on 
the interacting parameter e from the most recent data of 
SNe la (SNLS and Constitution samples), baryon acous- 
tic oscillation peak at z = 0.20 (2dFGRS) and z = 0.35 
(SDSS), and the Hubble evolution H{z). We have shown 
that negative and positive values for the interacting pa- 
rameter e are observationally allowed, which means that 
both an energy flow from dark energy to dark matter 
as well as a flow from dark matter to dark energy are 
possible. It is worth emphasizing that in our analysis 
the interacting parameter has been considered constant 
whereas in a more realistic case it must vary with red- 
shift. The theoretical and observational consequences of 
this more realistic interacting scenario, as well as a full 
comparison with the case discussed in the present analy- 
sis, will appear in a forthcoming communication [2l|. 
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